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Abstract. In this paper in space the problem of construction of optimal quadra- 

ture formula in the sense of Sard is considered and using S.L.Sobolev's method it is obtained new 
optimal quadrature formula of such type. For the optimal coefficients explicit formulas are obtained. 
Furthermore, the numerical results which confirm the theoretical results of this work is given. 
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1. Introduction. Statement of the problem 
We consider the following quadrature formula 

1 N 



j <p{x)dx=^Cp<p(xp) (1.1) 

0=° 



with the error functional 



N 

£(x)=e m (x)-Y,Cp5(x-x p ), (1.2) 

/3=0 

where Cp are the coefficients and xp are the nodes of formula (1.1), xp G [0,1], £[o,i](x) is the 
indicator of the interval [0,1], 5(x) is Dirac's delta-function, function tp(x) belongs to Hilbert space 
^'(O,!). The norm of functions in this space is defined by the following equality 



1/2 

rtxWt'^iO,!^ = I I {^ m \x) + ^ m ~ l \x)Ydx\ . (1.3) 




i 

2 



The difference 

N 



(£(x),v(x))= / (p(x)dx - C p ip(xp) = / £(x)<p{x)dx (1.4) 
o 0=° 

is called the error of the quadrature formula (1.1). The error of the formula (1.1) is linear functional 

in W^ m,m 1 ^*(0, 1), where W^™'" 1 (0, 1) is the conjugate space to the space W^™'" 1 (0, 1). 
By Cauchy- Schwartz inequality 

\(l(z),<p(x))\ < \W{x)\wt' m ~ l) (0,1)|| • U(x)\wt' m - lh (0,1)|| 
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the error (1.4) of formula (1.1) is estimated with the help of the norm 



i{x)\wt' m ~ 1] *(vx 



sup 

L(x)|w, ( " , ' m - 1) (o ) i)| 



\(e(x)M*))\ 



of the error functional (1.2). Consequently, estimation of the error of the quadrature formula (1.1) on 
functions of the space W^ 71 '™ ^(0, 1) is reduced to finding the norm of the error functional £{x) in 



the conjugate space W\ 



(m,m— 1)* 



(0,1). 



Obviously the norm of the error functional £{x) depends on the coefficients Cg and the nodes Xg. 
The problem of finding the minimum of the norm of the error functional £{x) by coefficients Cg and 
by nodes xg, is called the S.M.Nikolskii problem, and obtained formula is called optimal quadrature 
formula in the sense of Nikolskii. This problem was first considered by S.M.Nikolskii [15]. The problem 
were further investigated by many authors for various cases (see e.g. [2-5,16,33] and references therein). 
Minimization of the norm of the error functional £{x) by coefficients Cg when the nodes are fixed is 
called Sard's problem. And the obtained formula is called optimal quadrature formula in the sense of 
Sard. This problem was first investigated by A.Sard [17]. 

There are several methods for the construction of optimal quadrature formulas in the sense of Sard 
such as the spline method, the f(x)— function method (see e.g. [2], [24]) and Sobolev's method. Note 
the Sobolev method is based on the construction of a discrete analogue of a linear differential operator 
(see e.g. [27-29]). In different spaces based on these methods, the Sard problem was investigated by 
many authors (see, for example, [1,2,4,6-8,11-14,18-20,23-29,31,32] and references therein). 

Furthermore, explicit formulas for coefficients of optimal quadrature formulas for any m and for any 
number N + 1 of the nodes xg in the space L 2 "^ when the nodes are equally spaced were obtained in 

the works [11,18,19]. Here is the Sobolev space of functions, which m— th generalized derivative 
is square integrable. In the work [20] positiveness of the coefficients of optimal quadrature formulas 
in L 2 m ^(0, 1) space is investigated. 

It should be noted that in the works [8,12], using the ip— function method, investigated the problem 
of construction of optimal quadrature formulas in the sense of Sard which are exact for solutions of 
linear differential equations and several examples given for some number of the nodes. 

The main aim of the present paper is to solve the Sard problem in the space W^™ 1 '™ 1 ^(0, 1) using 
S.L. Sobolev's method for any number N + 1 of the nodes xg, i.e. finding the coefficients Cg satisfying 
the following equality 



£(x)\W, 



(m,m—l)* 



(o,i; 



inf 



£(x)|W 2 (m ' m - 1} *(0,l) 



(1.5) 



Thus, in order to construct the optimal quadrature formula in the sense of Sard in the space 



W. 



(m,m— 1) 



(0, 1) we need consequently to solve the following problems. 



Problem 1. Find the norm of the error functional £{x) of quadrature formulas (1.1) in the space 
W 2 (ro ' ro - 1) *(0,l). 

Problem 2. Find the coefficients Cg which satisfy equality (1.5) when the nodes xg are fixed. 

The structure of the present paper is: in the second section the extremal function which corresponds 
to the error functional £{x) is found and with its help representation of the norm of the error functional 
(1.2) is calculated, i.e. the problem 1 is solved; in the third section in order to find the minimum 
of the quantity \\£\\ 2 by coefficients Cg the system of linear equations is obtained for the coefficients 



of optimal quadrature formula in the sense of Sard in the space W. 



(m,m— 1) 



(0,1), moreover existence 
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and uniqueness of the solution of this system are proved; in the forth section explicit formulas for the 
coefficients of optimal quadrature formula of the form (1.1) are found, i.e. the problem 2 is solved; 
finally, in the fifth section the results of numerical experiments are given. 

2. The extremal function and the representation of the norm 
of the error functional £(x) 

In order to solve the problem 1, i.e. for calculation the norm of the error functional (1.2) in the 

space W^™ 1 '™ 1 ^*(0, 1), it is used the concept of extremal function of given functional. The function 
tpe(x) is called the extremal function for the functional £(x) (see, [28]), if the following equality is 
fulfilled 

{£{x),Mx)) = l^l^^-^^l)! • II^IW^'^O,!)!! • (2.1) 
This means that the extremal function ipe(x) gives the greatest error to the difference (1.4) in the 

TTr (m,m— 1) m -i \ 

space W$ (0, 1). 

Since the space W 2 (m ' m_1) (0,l) is Hilbert space then the extremal function tpe(x) in this space 
is found with the help of Riesz theorem about general form of a linear continuous functional on 
Hilbert spaces. Then for functional £(x) and for any ip{x) G W^" 1 '" 1 ^(0,1) there exists the function 
ipe(x) 6 W^ ri ' m ~ l \^, 1) for which the following equation is take placed 

(£(x),ip(x)) = <V>/(x), (2-2) 

where 

l 

y> t (x),<p{x)) = J (4 m) ( x ) + 4 m ~ 1) ( x )) (v (m) {x) + ^ m - l) (x))dx (2.3) 



is the inner product defined in the space 
Further, we will solve equation (2.2). 

o (oo) o (oo) 

Suppose (p(x) belongs to the space C (0, 1), where C (0, 1) is the space of functions, which are 
infinity differentiable and finite in the interval (0, 1). Then from (2.3), integrating by parts, we obtain 

l 

(Mx)M*)) = ("I)™/ (VF m) (z) -4 2m ~ 2 \x)) V>(x)dx. (2.4) 
o 

Keeping in mind (2.4) from (2.2) we get 

4 2m \x) - ^f m ~ 2) {x) = (-l) m £(x). (2.5) 

So, when <p(x) GC (0, 1) the extremal function ipe(x) is a solution of equation (2.5). But, we have 

[-IN O (OO) 

to find the solution of equation (2.2) when (p(x) G W^'™ 1 (0,1). Since the space C (0,1) is 
densely in the space W^ m ' m 1 - ) (0, 1), then we can approximate arbitrarily exact functions of the space 

( —1} ° ^ OC ^ (m —11 

yy(m,m J (o, 1) by a sequence of functions of the space (7 (0, 1). Next for any ip(x) G W^ 1 '" 1 (0,1) 
we consider the inner product (ipi(x),ip(x)) and, integrating by parts, we have 

m— 1 

M*)M*)) = E(- 1 )> (m " 1_s) (-) (^ m+s) (^) - 4 m+s ~ 2 \x)) \T=l+ 

s=l 
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1 

W m - l) (x) (4 m \*) + 4 m ~ 1) W>) l-o + ("I)™/ ^f m) (x) - ^ 2 \x))dx. 



Hence from arbitrariness <p(x) and uniqueness of the function ipn{x) (up to the function e~ x and 
polynomial of degree m — 2), taking into account (2.5), it must be fulfilled the following equations 

4 2m \x) - ^f m ' 2) (x) = (-l) m e(x), (2.6) 
^ { r +S \x)-4 m+s - 2 \x)) 1-1 = 0, (2.7) 

(^ m) («)+^ m - 1) ( a :))|^ = 0. (2.8) 

Thus, we conclude that the extremal function ipe(x) is the solution of the boundary value problem 
(2.6)-(2.8). 

The following holds 

Theorem 2.1. The solution of the boundary value problem (2.6)-(2.8) is the extremal function 
il>e(x) of the error functional i{x) and has the following form 

Mx) = (-l) m i(x) * G{x) + P m - 2 (x) + de~ x , 

where 

- g w^w) (2 ' 9) 

is a solution of the equation 

iP (2m \x) - V (2m " 2) (x) = S(x) (2.10) 

d is any real number, P m _2(x) is a polynomial of degree m — 2. 

Proof. It is known that general solution of nonhomogeneous differential equation consists on sum 
of a partial solution of nonhomogeneous differential equation and general solution of corresponding 
homogeneous differential equation. 

The homogeneous equation for differential equation (2.6) has the form 

^ m \x)-4 2m - 2 \x) = 0. (2.11) 
It is easy to show that general solution of equation (2.11) is 

P 2m _ 3 (x) + d 1 e x + d 2 e- x . (2.12) 
It is not difficult to verify that a partial solution of the differential equation (2.6) is 

(-i) m e( x )*G{x), 

where G(x) is a fundamental solution of equation (2.6) and is defined by (2.9) and is a solution of 
equation (2.10), * is operation of convolution, where convolution of two functions is defined as 

oo oo 

f(x)*g{x)= j f(x-y)g(y)dy = J f(y)g(x-y)dy. 



The rule of finding a fundamental solution of a linear differential operator 

d n d n 
dx n 1 dx r 



d n d"- 1 
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where a,j are constants, is given in [30, p. 88]. Using this rule, it is found the function G(x), which is 
a fundamental solution of the operator — f^ m -i an d have the form (2.9). 
Thus, we have the following general solution of the equation (2.6) 

M%) = (-1)"^) * G(x) + P 2m - 3 (x) + d ie x + <% x , (2.13) 
where P2m-z(x) = &2m-3^ 2m ~ 3 + ••• + b\x + bo is a polynomial of degree 2m — 3, d±, d 2 are constants. 

In order that in the space W^™'™ 1 ^(0, 1) the function tpi(x) will be unique (up to the function e~ x 
and polynomial of degree m — 2), it has to satisfies the conditions (2.7), (2.8). Here derivative is in 
generalized sense and 

^f\ x ) = (-l) m i(x) * G^ k \x) + p£_ 3 (z) + d x e x + {-l) k d 2 e-\ 



where 



< e x - e~ x 



G(fc)(s) = signx 



k = 1,2, ...,2m - 1, 

j,2m- 3— k 

~ (2m-3-fc)! ~ 

' x 2 



3! 



— x when k is even, 



g x _|_ g — x „2m—3—k 



— — — ... — 1 when k is odd. 

(2m- 3- jfe)! 2! 



From conditions (2.7) for s = m — 1, taking into account 

e x + e~ x 



i)f m - 1] {x) = (-l) m £(x) * jsignx 



+ d±e x - die x , 



(2m-3), 



-l) m £(x) * 



signr 



e x + e x 



- 1 



+ (2m - 3)!6 2m _ 3 + d x e x - d x e 



we get 



4>f m - l \x) - 4 2m ~ 3 \x) = (-i) m £(x) * 



signs 



- (2m - 3)!&2m-3 = 



(- ' )'" (%), Sign( o ) ■ (2m- Win,-*- 



Hence for x = 0: 



4 2m - l) (o)-c"'" o; (o) = (-ir(%), 



, (2m-3) 



sign(-y) 



- (2m-3)!6 2 m- 



(-iy 



(£(y),l)-(2m-3)!62m-3 = 0, 



for x = 1: 



(2m-l), 



(l)-4 2m - 3) (l) = (-ir(^(y) 



sign(l - y) 



- (2m - 3)!6 2m - 3 = 



("l) r 



1) - (2m - 3)!62m-3 = 0. 



Then 



6 2m _ 3 = and (€(s/), 1) = 
When s = m — 2 from conditions (2.7), taking into account 

,(2m-2) 



(2.14) 



'( x ) = (-l) m ^(x) * <^ signs • 



+ die x + die x , 
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^-4) (I) = ( _ ir<N(x) , |2|£ [£Ll£I - i J + 
+ (2m - 3)!6 2m -3^ + (2m - 4)!6 2m _ 4 + die 21 + d 2 e~ x , 

we have 

Vf m " 2) (x) - V? m ~ 4) (*) = (-l) m ^(x) * - (2m - 3)!6 2m _ 3 x - (2m - 4)!6 2m _ 4 . 

Hence, taking account of (2.14), for x = 0: 

("l) m ^sign(0 - y)\ - (2m - 4)!6 2m _ 4 = 



(-1)' 



(%),y)-(2m-4)!6 2m _ 4 = 0, 



for x = 1: 



(-])'" (£(y),i_l s ign(l-y)) - (2m - l)!ft 2 ,„. , 



("l) f 



(£(y),y)-(2m-4)!6 2m _ 4 = 0. 



Hence 



&2m-4 = and (£(y),y) = 0, 
and so on, continuing by this manner, for s = m — 3, m — 4, 2, 1 we obtain 

6 m - 2+s = 0, Wy),j/ ro - 1 -) = 0. 

Combining (2.14), (2.15), (2.16) we get 

&2m-3-s = 0, s = 0,1,..., m-2, 

(£(y),y s ), s = 0,1,..., m-2. 
From the condition (2.8), keeping in mind (2.17), we have 

^ m) (x) = {-l) m l{x) * G (m) (x) + d ie x + (-l) m d 2 e- x , 
^'^{x) = {-l) m l{x) * G^-^ix) + d ie x + {-l) m ~ 1 d 2 e 

Hence 



(-l) m £(x) * 



signx 



x 



m-2 



„m-3 



For x = 0, taking into account (2.18), 



^ m) (o) + vT~^(o) = (-i) m (%)- 



( ,(m-l)/ n x _ 



(m-2)! (m-3)! 

sign(0 - y) 



(-v) 



m— 3 



(m - 3)! 



2 

_ l) m 



x-ljj+d ie x . 
<-y) m - 2 



-v 



(m-2)! 
{£ N (y),e-y)+d 1 =0. 



For x = 1: 



^(i) + ^(i) = (-i) m U(y), 



,(m-l) 



sign(l - y) _ / _ (1 - y) 



m-2 



(m - 2)! 



(2.15) 
(2.16) 

(2.17) 
(2.18) 



(i-y) 



m— 3 
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(-l) m e 



(m-3)! 



(l-y)-l +d ie 



(€jv(y),e- y ) +die = 0. 



Whence 

di = 0, (2.19) 

(£(y),e-y)=0. (2.20) 

Taking into account equalities (2.17)- (2.20) and denoting d.2 = d, we get the statement of the theorem. 
Theorem 2.1 is proved. □ 

The equalities (2.18) and (2.20) mean that our quadrature formula is exact for the function e~ x 
and for any polynomial of degree up to m — 2. 

Now, using the result of theorem 2.1, we immediately obtain the representation of square of the 
norm of the error functional (1.2) 

" N N 
. 0=0 7=0 

1 



e(x)\wt' m ~ 1] *(o,i) = (£(x),Mx)) = (-i) 



AT 1 11 

-2y^ C/3 J G(x-xp)dx + j j G(x-y)dxdy 



(2.21) 



/ ?=0 o o o 

Thus the problem 1 is solved. 

Further in sections 3 and 4 we solve the problem 2. 

3. The system for the coefficients of optimal quadrature formulas in the space 



W. 



(m,m—l) 



(0,1) 



Assume that the nodes xp of the quadrature formula (1.1) are fixed. The error functional (1.2) 
satisfies conditions (2.18) and (2.20). The norm of the error functional i{x) is multidimensional 
function with respect to the coefficients Cp {(3 = 0,N). For finding the point of conditional minimum 
of the expression (2.21) under the conditions (2.18) and (2.20) we apply the method of undetermined 
multipliers of Lagrange. 

We denote C = (Co, C±, Cn) and A = (Ao, Ai, A m _i). 

Consider the function 

ro-2 

*(C, A) = \\l{x)\\ 2 - 2(-l) m £ K (£(x), x a ) - 2(-l) m A ro _i {£(x), e~ x ) . 

a=0 



Equating to partial derivatives of ^(C, A) by Cp {(3 = 0, N) and Ao, Ai, A m _i, we get the following 
linear system 

N m-2 

£ C^G{xp - x 7 ) + J2 + x rn-ie~ x e = f m (xp), P = OTW, 

7=0 a=0 

N 1 

V C 7 x Q = — — , a = 0, 1, m - 2, 
„ a + 1 



7=0 



N 



^C 7 e-^ = l-e-\ 

7=0 



(3.1) 
(3.2) 
(3.3) 



8 KH.M.SHADIMETOV, A.R.HAYOTOV 

where G(x) is defined by equality (2.9), 



fm(xp) = / G{x - xp)dx. (3.4) 

J 



The system (3.1)-(3.3) has a unique solution and this solution gives minimum to ||^(x)|| 2 under the 
conditions (3.2), (3.3). Uniqueness of the solution of such type of systems is discussed in [28, 29]. 

It should be noted that existence and uniqueness of optimal quadrature formulas in the sense of 
Sard is also investigated in [12]. 

Now in (2.21) we will do change of variables Cp = Cp + C\p. Then (2.21) and the system (3.1)-(3.3) 
have the following form: 



(-1)' 



N N N 



£ CpC^Gixp - x 7 ) - 2 ( c /3 + c i,p) / G(x - xp)dx+ 



7=0 0=0 

1 1 



N N 

+ 

/3=0 7 =0 £ o 



^ (2C^C li7 + Ci ija Ci, 7 ) G(x p - x y ) + J J G(x - y)dxdy 



(3.5) 



N m-2 

YC^G{xp-x^ + Y, A a x^ + A m _ ie -^ =F m (xp), /3 = 0^V, (3.6) 

7=0 a=0 



N 



^C 7 x^ = 0, a = 0,m-2, (3.7) 

7=0 

TV 

E^ = 0, (3-8) 

7=0 

N 

where F m (xp) = fmixp) — C\ n G{xp — x 7 ), C\$ is a partial solution of equations (3.2), (3.3). 

7=0 

Hence we directly get that the minimization of (2.21) under the conditions (3.2), (3.3) with respect 
to C/3 is equivalent to the minimization of expression (3.5) with respect to Cp under the conditions 

(3.7) , (3.8). Therefore it is sufficient to prove that the system (3.6)-(3.8) has a unique solution with 
respect to unknowns C = (Co, Ci, C/v), A = (Ao, Ai, A m _i) and this solution gives conditional 
minimum to the expression ||^|| 2 . 

From the theory of conditional extremum it is known the sufficient condition in which the solution 
of the system (3.6) - (3.8) gives conditional minimum to the expression ||£|| 2 on the manifold (3.7), 

(3.8) . It consists on positiveness of the quadratic form 

N N 



*( c ) = EE^r^ (3.9) 

/3=0 7 =0 ^P^l 



on the set of the vectors C = (Co, C\, Cat) under the condition 

SC = 0, (3.10) 
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where S is the following matrix of equations (3.7), (3.8): 

/ 1 1 ■■ 



x 



Xl 



„m—2 vn— 2 



"0 



X 



-Xo 



-Xl 



1 \ 



X N 



X 



m-2 
N 
-x N 



J 



We will show that in our case this condition is fulfilled. 

Theorem 3.1. For any nonzero vector C G R N+1 lying in the subspace SC = 0, the function 3>(C) 
is strictly positive. 

Proof. Using the definition of the function \I/(C,A) and equations (3.5), (3.7), (3.8) from (3.9) we 
get 



N N 



d> (C) = 2(-l) m £ £ G(^ - x 7 )^C 7 . 
/3=0 7=0 

Consider the linear combination of delta functions 



N 



5 u (x) = V2Y,C{,6(x-x f) ). 



(3.11) 



(3.12) 



By virtue of the condition (3.10) this functional belongs to the space W 2 (m ' m - 1) *(0,l). So, it has 
the extremal function u-^{x) G W 2 ^ m ' m '~ 1 ' ) (0, 1) which is a solution of the equation 



^-^)« c w = (-ir%w. 

As U(=t(x) we can take linear combination of shifts of the fundamental solution G(x): 

N 

uc(x) = ^(-l) m ^^G(x-x^). 
The square of its norm in the space W^ m ' m-1 ^(0, 1) coincide with $ (C), i.e. 

2 N N 

c(x)\wt m - l) (0,l) = (%(x),n n (x)) =2(-ir££^C f 7 G(x / 3-x 7 ). 



(3.13) 



=0 7=0 



Hence clearly that for nonzero C the function <3? (C) is strictly positive. 
Theorem 3.1 is proved. 



□ 



If the nodes xq,xi, ...,xn are selected such that the matrix S has right inverse matrix, then the 
system (3.6) - (3.8) has a unique solution. Then the system (3.1)-(3.3) also has a unique solution. 

Theorem 3.2. // the matrix S has right inverse matrix, then the main matrix Q of the system 
(3.6) - (3.8) is nonsingular. 

Proof. We denote by M the matrix of quadratic form =■ 2 > $ (Cj, where <I> (Cj is defined by 
equality (3.11). It is known that if homogenous system of linear equations has only trivial solution 
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then corresponding nonhomogeneous system has a unique solution. Consider homogeneous system 
corresponding to the system (3.6) -(3.8) in the following matrix form: 



Q 



c 

A 



M S* 
S 



c 

A 



0, 



(3.14) 



where S* is the transposed matrix to the matrix S. 

We verify, that unique solution of (3.14) is identical zero. Suppose C, A is the solution of (3.14). 
Consider the function <%(#), which is determined by equality (3.12). As the extremal function for the 
function <%(x) we take the following function: 



N 



m—2 



u-, 



T ( x ) = v^-ir E c p g ( x - x ^ + E A « x ° + A ™- ie " 



a=0 



This is possible because the function u-^(x) belongs to the space is a solution of 

equation (3.13). First A + 1 equations of the system (3.14) mean that Uq(x) takes values at all 



nodes xp. Then for the norm of the functional <%(x) in W 2 



(m,m— 1)* 



(0, 1) we have 



N 



Sc(x)\wt' m - 1] *(0,l) = (%W,«cW) =2(-l) m Y l Cpu c {xp) = (3.15) 



P=0 



On the other hand taking into account equations (3.7), (3.8) we get 

N N 

r {m,m—i)* 



6c(x)\W? 



"(0,1) 



/3=0 7=0 



(3.16) 



From (3.16) we conclude that (3.15) is possible if C = 0. Then from first + 1 equations of the 
system (3.14) we obtain 

S*X = 0. (3.17) 
By assertion of the theorem, the matrix S has right inverse matrix, then S* has left inverse matrix. 
Then from (3.17) we conclude that the solution A also is equal to zero. 

Theorem 3.2 is proved. □ 



From (2.21) and theorems 3.1, 3.2 it follows that in fixed values of the nodes xp the square of 
the norm of the error functional £(x) being quadratic functions of the coefficients Cp has a unique 

o 

minimum in some concrete value Cp = Cp. 

o 

As said in the first section the quadrature formulas with the coefficients Cp (ft = 0, A), corre- 
sponding to this minimum in fixed nodes xp is called optimal quadrature formula in the sense of Sard 

o 

and Cp (ft = 0, A) are called optimal coefficients. 

o 

Below for convenience the optimal coefficients Cp remain as Cp. 

4. Coefficients of optimal quadrature formula in the sense of Sard 

In the present section we will solve the system (3.1)-(3.3) and will find explicit formula for the coef- 
ficients Cp. Here we will use similar method suggested by S.L.Sobolev [27] for finding the coefficients 

of optimal quadrature formulas in the space L^^O, 1). Here mainly is used the concept of functions of 
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discrete argument and operations on them. Theory of discrete argument functions is given in [28,29]. 
For completeness we give some definitions about functions of discrete argument. 

Assume that the nodes xp are equal spaced, i.e. xp = h/3, h = jj, N = 1, 2, .... 

Suppose that <p(x) and ip(x) are real-valued functions of real variable and are defined in real line R. 

Definition 4.1. A function <p(h(3) is called function of discrete argument if it is given on some set 
of integer values of f3. 

Definition 4.2. By inner product of two discrete functions <p(h/3) and i/>(h/3) is called the number 

oo 
/3=-oo 

if the series on the right hand side of the last equality converges absolutely. 

Definition 4.3. By convolution of two discrete functions (f(h(3) and ip(h/3) is called the inner 
product 

oo 

l p(h/3)*4,(hp) = [<p(h 1 ),4>(hP-h 1 )}= J2 vW-Mhp-h-y). 

7=— oo 

Now we turn on to our problem. 

Suppose that (7/3 = when (5 < and f3 > N. Using above mentioned definitions the system 
(3.1)-(3.3) we rewrite in the convolution form 

G{hf3)*Cp + P m -2{hf3) + de- h P = f m (hl3), (3 = 0, 1, N (4.1) 

C p = 0, when (3 < and /3 > N, (4.2) 

N 1 

VC r (^r = — , a = 0,l,...,m-2, (4.3) 
ft a + 1 

N 



^Cp-e-W = l-e~\ (4.4) 



where 



. (hR . e h ? + e-^ + e 1 -^ + e^" 1 - 4 ^ (^/3) 2fc + (1 - fe/3) 2fc 
/m^J- 4 ^ 2 (2A;)! • (4.5) 

Consider the following problem. 

Problem A. Find the discrete function Cp, polynomial P m -2{hj3) of degree m — 2 and unknown 
constant d which satisfy the system (4-l)-(4-4) f or given f m (h/3). 

Further we investigate the problem A and instead of Cp we introduce the functions 

v(h/3) = G(h/3) * Cp (4.6) 

and 

u(hP) = v {hp) + P m - 2 {hp) + d e~ h(3 . (4.7) 

In such statement it is necessary to express the coefficients Cp by the function u(h/3). For this we 
have to construct such operator D m (hf3) which satisfies the equality 

D m (h0)*G(hP)=6(hp), (4.8) 

where 5 (h/3) is equal to when (3^0 and is equal to 1 when f3 = 0, i.e. S(h/3) is the discrete 
delta-function. 
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In connection with this in [21,22] the discrete analogue D m (hf3) of the operator — d Ji m -i ■, 
which satisfies equation (4.8) is constructed and its some properties are investigated. 
The following theorems are proved in the works [21,22]. 

Theorem 4.1. The discrete analogue of the differential operator dx - lm — dx ^ m -i satisfying the 
equation (4-8) has the form 



where 



( m-l 

E A k \f~\ 

k=l 



D m (h(3) 



m—l 



(2m-2) 
^2m-2 



-2e h + Z A k , 
k=l 
m—l 

k=i 



\P\ > 2, 
l/3| = 1, 
P = 0, 



C = 1 + (2m - 2)e n + e M + 



h . (2m-2) 
2h , e P2m-3 



_(2m-2) 
Plm-2 



2m- 2 



2(1 - A fe ) 2 "- 2 [A fe (e 2 ^ + 1) - e>*(A 2 + l)bfc 2) 



^fe - P 2m-2( A fe) 

^2m- 2 (A) = ^ pf m - 2) A s = (1 - e 2h ){l - \) 2m - 2 - 2(\(e 2h + 1) - e h {\ 2 + l))x 



s=0 



fr 3 (l - A) 2m ~ 



•6 



_ A) 2m ~ 4 + ^iZ_v E 2 (X) + ... + " ^ 2m ~ 41 



/ l 2m - 3 ^ 2m _ 4 (A) 



(4.9) 



(4.10) 
(4.11) 



3! (2m -3)! (412) 

V2m-2' '> are the coefficients of the polynomial V 2m - 2 (\) defined by equality (4.12), A& are 

i/ie roots o/ £/te polynomial 7 ? 2m-2(A) which absolute values less than 1, Ek(X) is the Euler-Frobenius 
polynomial of degree k (see [29]). 



Theorem 4.2. T/ie discrete analogue D m (h(3) of the differential operator dx - 2m — dx 'i m ~i satisfies 
the following equalities 

1) D m (hP)*e h f > = 0, 

2) D m {h(5) * e-W = 0, 

3) D m {hp) * (h/3) n = 0, n < 2m - 3, 

4) D m {hp)*G{hp) = S{hp), 

here G{hj3) is the function of discrete argument corresponding to the function G{x) defined by equality 
(2.9) and 5{hj3) is the discrete delta function. 



Then taking into account (4.7), (4.8) and theorems 4.1, 4.2, for optimal coefficients we have 

C(3 = D m (h(3)*u(hf3). (4.13) 

Thus if we will find the function u{hf3) then the optimal coefficients will be found from equality 
(4.13). 

In order to calculate the convolution (4.13) it is required to find the representation of the function 
u{hf3) for all integer values of /3. From equality (4.1) we get that u(hf3) = f m {h[3) when h/3 € [0, 1]. 
Now we need to find the representation of the function u(h/3) when j3 < and (3 > N. 
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Since Cp = when h/3 ^ [0, 1] then 

Cp = D m (hf3) * u(hp) = 0, hfit [0, 1]. 

Now we calculate the convolution v{h[i) = G{hj3) * Cp when hp £ [0, 1]. 

Suppose < then taking into account equalities (2.9), (4.2), (4.3), (4.4), we have 

oo 

v(hp) = G(hp) *Cp= c i G ( h P ~ h i) = 

7=— oo 

_ V-„ signW - h '<) ( ehf, ~ hl - e- h f+ h ~< "^j (h/3 - /i 7 ) 2t -' \ 

<= w ^ „ , e-w " „ h i" ^ (hp - At) 2 *- 1 



7=0 7=0 7=0 fc=l 

^ e-^A fc 1 N 



= -^-(i-e- 1 ) + ^ T -E^ 7 + ^E^ 

7=0 7=0 



' 1 (fe/3) 2fc - 1 -"(-/l 7 ) a 



V 



(2k- 1 - a)! • a! 

fc=l Q=0 V ^ 



™ 1 ^ 2 (/ t/ 3) 2fc ~ 1 ~ Q (-/ t7 )" ^ 2 ^ (/i/3) 2fc - 1 ^(-/i 7 ) a 
+ 2^ 2^ (2fc - l-a)!-a! + ^ ^ (2fc - 1 - a)! • a! 

fe= ^ m + l j a= Q V ' /;= ^ m + l j g= m -l y ' 

p hfi p -hj3 N i m ~ 1 2k ~ l ( h o\2k-l-a( i\« N 

-V^+VE^+I E E fb^E^r+ 

7=0 fc = j22±ij a=m— 1 ' 7=0 



[ m + l 

i r ^ (hp)' 2k -^ a (-i) a (hpy 2 *- L - a (-i) a 



-til —1 \ 

+ 2 2_> 2_> (2fe-l- a )!.(a + l)! + ^ ^ (2fc - 1 - a)! • (a + 1)! 

\ i n 1 fe= ^ m+l j a=0 v I \ ' J j 



k=l a=Q 



Thus we get 



where 



= -^-(1 - e" 1 ) + De"^ + Q( 2m ~ 3 )(/i/3) + Q m - 2 {hp). 

p hj3 

v (hp) = - — (1 - e^ 1 ) + De-W + Q( 2m - 3 )(/t/3) + Q m - 2 {hp). (4.14) 

([m±ll_l 
& ^(2*-l-a)!-(a + l)! + 

+ L L (2fe _ 1 _ a)! . (a + 1)! I ( 4 - 15 ) 

fc= [Z2+i] a=0 V 
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is the polynomial of degree 2m — 3 with respect to (h/3), 



Q~*W-\ E E (4.16) 

fc=^Hi±I] a=m— 1 v ' 7=0 

is unknown polynomial of degree m — 2 with respect to (/i/3), 

1 N 

7=0 

Similarly, in the case f3 > N for the convolution v (h/3) = G(h(3) * Cp we obtain 

v(h/3) = e —(l - e- 1 ) - De-W - Q^ 2m - 3 \h(3) - Q m - 2 (hf3). (4.18) 

We denote 

QmUW) = Pm-2(h/3) + Q m -2(h(3), aT = d + D , (4.19) 

QmUhP) = Pm-2(hP) - Q m -2(h(3), a + = d - D . (4.20) 
and taking into account (4.14), (4.18), (4.7) we get the following problem. 
Problem B. Find the solution of the equation 

D m (h/3) * u(h/3) =0, hpt [0, 1] (4.21) 

having the form: 

-4^(1 - e" 1 ) + a-e-^ + Q( 2m ~ 3 ) (hp) + Q { ~l 2 (hp), P < 0; 

/m(^), < /? < N; (4.22) 



u(/i/3) 



^(1 - e" 1 ) + o+e-^ - Q( 2m " 3 )(/ l /3) + q£2 2 (/*0), /3 > JV. 



.Here Q[ n }_ 2 (hf3) and Q^}_ 2 (h/3) are unknown polynomials of degree m — 2 mi/t respect to h/3, a and 
a + are unknown constants. 

If we find Q { ~± 2 (hp), Q { +l 2 (h(3), a~ and a+ then from (4.19), (4.20) we have 
P m - 2 (hP) = \ (Q ( m l 2 (h{3) + Q ( +l 2 (hp)) , d = \(a~ + a + ), 

Qm-2(h(3) = \ (Q { m l 2 (h[3) - Ql%(hp)) , D = l -(a- - a+), 

Unknowns Q^2 2 (/i/3), Q^}_ 2 (h/3), a~ and a + can be found from the equation (4.21), using the 
function D m (h/3). Then we can obtain explicit form of the function u(hf3) and find the optimal 
coefficients Cp. Thus, the problem B and respectively the problem A can be solved. 

But here we will not find Q^ m - 2 (hP), Q^l 2 (hp), a~ and a + . Instead of them, using D m (h(3) 
and u(h/3), taking into account (4.13), we will find the expressions for optimal coefficients Cp when 
p = l,...,N-l. 

We denote 

ak = J^Yl X k {- V (1 " ^ + Q (2m ~ 3) (-hl) + QtU-hl) + a~e^ - f m (-h 7 )j , (4.23) 



OPTIMAL QUADRATURE FORMULAS 15 
a oo / ^7+1 \ 

(4.24) 

Here A& are the roots and p is the leading coefficient of the polynomial T > 2m-2(X) of degree 2m — 2 
defined by equality (4.12) and |A&| < 1. The series in the notations (4.23), (4.24) are convergent. 
The following is true 

Theorem 4.3 (Theorem 3 of [23]). The coefficients of optimal quadrature formulas in the sense of 
Sard of the form (1.1) in the space W^ 71 '™ ^(0, 1) have the following form 



m—l 



Cp = D m (h(3)*f m (h(3) + Y,(ak4 + b k \%- fS ), = 1, 2, N - 1, (4.25) 

k=l 

where and b^ are unknowns and have the form (4-23) and (4-24) respectively, A& are the roots of 
the polynomial 7- > 2m-2(A) which defined by equality (4-12) and \\k\ < 1. 

Proof. Suppose (3 = 1,N — 1. Then from (4.13), using (4.9), (4.22) and theorem 4.2, we have 

oo 

Cp = D m (h{3)*u(h(3)= D m {hp-h 1 )u{h 1 ) = 

7=— oo 

— 1 N oo 

= D m (hP-h 7 )u(h 7 )+J2 D m(hP-h 1 )u(h 1 )+ £ Dm(hP - h^uihj). 

7=— oo 7=0 7=W+1 

Hence using the definition of convolution of discrete functions we get 

Op = D m (hp) * f m (hP) + 

oo m—l A / 



oo ill j. . , 

E E ^ A f +7 (" V (1 " 6_1) + ^ (2m ~ 3) ("^) + QnU-hl) + a-e^ - / m (-fry)J + 

oo m—l . / 

+ EE - e" 1 ) - Q( 2 - 3 )(l + fry) + QW(1 + fry) + a+e- 1 "^- 

7=1 fc=i Afc P V 4 



-/m(l + /»7) I- 

Whence taking into account the notations (4.23), (4.24), we get (4.25). 

Theorem 4.3 is proved. □ 



From theorem 4.3 it is clear that in order to obtain explicit forms of the optimal coefficients Cp 

in the space W^ rt " m ^(0, 1) it is sufficient to find au and bk (k = l,m — 1). But here we will not to 
calculate series (4.23) and (4.24). Instead of them substituting the equality (4.25) into (4.1) we obtain 
identity with respect to (h/3). Whence, equating corresponding coefficients the left and the right hand 
sides of equation (4.1) we will find a& and bk- And the coefficient Co and Cn will be found from (4.3) 
when a = and (4.4), respectively. Below we will do it. 

It should be noted that the cases m = 1 and m = 2 are solved in the work [23] and the following 
theorems are proved. 
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Theorem 4.4 (Theorem 4 of [23]). The coefficients of optimal quadrature formulas of the form 
(1.1) with equal spaced nodes in the space W^'^iO, 1) are expressed by formulas 



C ft 



e h +l ' 
2(e h -l) 



P = 0,N, 



p = l,N-l, 



where h= l/N, N = 1,2,.... 

Theorem 4.5 (Theorem 5 of [23]). The coefficients of optimal quadrature formulas of the form 

(2 1) 

(1.1) with equal spaced nodes in the space W 2 ' (0, 1) are expressed by formulas 

' l--h--K(h)(\ l -\f), (3 = 0, 



h + K(h) ((e h - Ai)A? + (1 - Aie h )Af- /3 ) , = 1, N - 1, 



where 



Ai 



/i(e 2/l + 1) - e 2/l + 1 



(2e fe -2-/ie fe -/i)(Ai-l) 



JV+In 



2(e> l -l) 2 (A 1 + Af 
(e' 1 - l)V/i 2 (e ft + l) 2 + 2h(l-e h ] 



l-e 2h + 2he h 



, |Ai|<l, 



h = l/N, N = 2,3,.... 

The main goal of the present section is to solve the system (4.1)-(4.4) for any m > 2 and any natural 
N, N > m. As mentioned above for this sufficiently to find and (k = l,m— 1) in (4.25). 
The main result of the present paper is the following theorem. 

Theorem 4.6. XTie coefficients of optimal quadrature formulas of the form (1.1) with the error 
functional (1.2) and with equal spaced nodes in the space W^ 71 '™ ^(0, 1) when m > 2 are expressed 
by formulas 



C 



-i-h 



m—l 



t+ E 

k=l 



\ k (e h -e)+\l(e-l)+\% +1 (l-e h ) | A£ r+1 ( e h - e )+Aff(e-l)+A fc (l-e h ) 



(e-l)(l-A fc )(e h -A fc ) 



(e-l)(A fc -l)(A fc e fe -l) 



m—l /a o\ 

C/j = fc + £ (afcAf + /3 = l,iV-l, 

fc=l v 7 



m—l 



n e h h+l-e h | 



afc- 



A fc (e-e h + 1 )+Aj r (c h + 1 -e h )+Aj r+1 ( C h -e) 



fc=l 



(e-l)(l-A fc )(e h -A fc ) 



.JV+l 



(e-e' I + 1 )+A 2 (e h+1 -e' l )+A fc (e* 1 -e) 
(e-l)(l-A fe )(l-A fce h) 



where and (k = l,m — 1) are defined by the following system of 2m — 2 linear equations 



Afc 



m—l 

afc (A fc -l)(A fc -< ; i 
m—l \N+l 



m—l 

k-) + £ h 

k=l 
m—l 



h-2 



(A fc -l)(A fc e"-l) _ 2(e"-l) + (e*-l) a ' 



A "" *" A fe 

fc £ a k (Afc _ 1 ) (Afc _ eh) + 6fe (A fc -l)(A fc e*-l) 



2- 



m—l 

E a fe 
k=i 

m—l 

E 

fc=i 



J fe 2i-2 2 ' 2 AfcA 'Q 2 '- 2 

V- fe 2 '- 1 ^ A^A^O 2 '- 1 

i=i (2 ' _1)! ,4i 



m—l 

+ E h 
k=l 

m—l 

+ E h 
k=i 



3 

E 



h-2 i _ 
2(e h -l) (e"-l) 

(21^27 2- (1-A fc )* +1 



h 2i-2 Ai v+ 'A l 2; - 2 



E 



h 2i-l 



2i— 2 ,jv+i 

E ^ 
i=i 

2j-l A JV+ lAlQ 2!-l 



i (21-1)1 g (l-A fc )'+i 



0, J = 2,[f]; 



E 

Z=l 



fe 2 'B 
(20 



21 „■ 
I J J 



m— 1 1 
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m— 1 

fc=i 



f fr C * f ^ + ' A ' ' _ fr f ^L Al ° J 



i=i t=i i=i 



m—l 

k=i 



1=1 



i=l 



^ (Afc-l)^ 1 
1=1 

.Here are the roots of the polynomial (4-12) and \Xk\ < 1, Bj are Bernoulli numbers. 
In the proof of theorem 4.6 we use the following formulas from [10]: 

n—l -, k / \ i „ k 



7=0 H i=0 v y/ y i=0 v y/ 



/ l7=nv 



(4.26) 



where A*0 fe = Yl C\l k , A l ^ k is the finite difference of order i of j k , 

i=i 

and from [9]: 

13-1 k+1 7i r> 



§iK* + i-j)J 



where B^+i-j are Bernoulli numbers, 



A a a^ = E ( p ) A Q 0%^. 



p=0 

Proof of theorem 4.6. For convenience we denote 

T = D m (hp) * f m (hp). 

Now we consider equality (4.1) 



N 



E c , 

7=0 



sign(/i/3 - fry) / e ft P 



,hP-h-y _ Ji-y-hp 



m—l 

E 

fe=i 



fog ~ M 
(2fc-l)! 



2fc-l 



(4.27) 



(4.28) 



(4.29) 



where /3 = 0,1,..., TV. 
We denote 



+ P m _ 2 (fr3)+de-^ = / m (fr£0, 

(4.30) 



TV 



<7(fr3) = ^C 7 

7=0 



sign(fr3 - fry) / e^~ fe T - e fe T~^ _ y! (fr3 - fry) 2 ^ 1 ' 
9 9 2^ 



fc=l 



(2&-1)! 



(4.31) 



The expresson g{hj3) we rewrite in the following form 



JiB 



9{hf3)= Co 



E 



'hd-hf _ e hj-h/3 m ~ 1 



7=1 



E 

fc=l 



(2fc - 1)! 



(fr/3 - fry) 2 ^ 1 
^ (2fc- 1)! 



7=0 



fe=l 



(4.32) 
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Further we denote 



13-1 

g 1 (hf3) = Y / C 1 



7=0 



fc=l 



(2k- 1)! 



52 M 



1 



2 

7=0 



k=l 



(2fe-l)! 



e ^7 _ e -^7 



-E 



Firstly we consider g\(hj5) and rewrite (4.33) in the following form 

7=1 

Hence using (4.25), (4.29) we get 

/8-1 / m-l 



1 (2fc-l)iy 



p-i / m-l \ /j 7 



7=1 



fc=l 



0-1 



m—1 



7=1 
m-lP-1 



^-/3+7 



fc=l 
m—1 



-/17 



EE 



£=1 7=1 



fc=l 



2 

(fry) 2 *" 1 
(2£-l)!" 



From here taking into account (4.26), (4.27) and after some simplifications we have 



9i (h/3) 



T(e h - e h P) 



m—1 



—mv + E h 



Af e h - \ k eW 



2{l-e h ) 



e" - Af 



fc=i 



2{e h - 1) ^ 



afc 



2(A fe -e fc ) 2{l-\ k e h ) 
Aj-Afce*"^ . A^ +1 -A{V^' 



E 



m—1 



E 



€=1 
_2 m—1 



fc=l 

Thr\hp) 

~ w 



2(X k e h - 1) 



2(e^-A fc ) 



l=i j 



E 



,JV-fl 2^-1 



A 



,3+1 2^-1 A i Q 2£-l 



1-A fc 



E i 



i=0 



a k 



Afc 



Afc 



A fe 2 ^ A*/^" 1 



1)* A fc 



/- 

(A fc 
i=o v ft 



I) 1 



+ 



Afc 



A * 2£-l 



1-A fc 



2^-1 



Ei 



i=0 



Afc 



Afc 



AJ/3 2£-l 



(4.33) 
(4.34) 



(4.35) 



In the last expression of gi(hfi) the coefficients of A^ and 13 are the values of the polynomial 
^ , 2m-2(A) which defined by equality (4.12) at the Afc. Since Afc are the roots of the polynomial (4.12), 
then the coefficients of A^ and A^ _/3 are zero. Then from (4.35) for gi(h(3) we get 

m—1 



9iW) = — 



T Afc A^ \ 
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-h/3 



Te 



h m ~ l 



k=l 



N „h 



Al v e 



+ 



2£-l 



+ 



'VV, Uhi)^ T^,2*-i^F B n _ 3 



2{l-e h ) 2(1 -e 



E 



E E 



^■!(2^-i) 



m— 1 



m—l 



A fe 2 ^ A^ 2 ^ 1 . ^ ft*" 1 



m—l 



,N 2^-1 



+ E (2^-1)! E^ Afc _i E (A fc -1)< + E (2^-1)! ES-A fc E (x 

<f=l v y fc=l K i=0 v ; 1=1 K ' k=l i=0 v 

Finally taking into account equality (4.28) we have 



Afe 



A l /3 



ffl(^) 



rp m-1 

oh _ i + E 

fc=l 



Te h 



m—l 



Ji 



e h X k 



Te h 



X k -e h 

21-1 



2(1 - e h ) 2(1 - e* 



k=i 



+ 

B 2 £-j 



f3 j + 



m— 1 



A 



— 1 A^O 2 ^ -1 " 



+ E (2 ^ _ i)i E cit-ifi 3 E °* Afc _ i E (A _ 

£=1 v ; j=0 k=l K i=0 y K ' 

m—l 7 2£— 1 2<?— 1 m—l * jy 2^—1 ✓ < \ j 



+ E^E^E^E(^ 

Now using the binomial formula and equalities (4.3), (4.4) from (4.34) we obtain 



92{W) = - 



1 / e h P 



N 



7=0 



1 

~ ~2 

m—l m— 2 

+ E E 



2 

e hf3 



Ec, 



AT m-1 



E<V--E^E 



(^ - h 7 ) 



2m- 1 



7=0 7=0 fe=l 

r m+l 1 -i 

' L— a— J— 1 2fc— 1 /Lm2 fc-l-a 



2 ^-" 7 

7=0 



E E 



{2k- 1)! 



(/i/3) 2fe 



2fc-l-Q; 



■1)' 



+ E E 



(2fe-l-a)!(a + l)! 



\ fc=l a=0 
m-1 2fe-l „so fc _i_„ / AT 



^ (2/c - 1 - a)\(a + 1)! ^ (2k - 1 - a)\a\ ^ 

k= ™±1 a=0 v l\ J k _ m+i Q = m -1 v ; 7=0 

Using the binomial formula in equality (4.5) for f m {h(3) we have 



Ec^r 



m—l 



i [ ^~ 2 ^ m 2k - > 



(4.36) 



(4.37) 



t (2*)! 



(-1) Q 



i ^ ^ 2 (^^-^"(-l)" i ^ '^-r (^^-^"(-i)" 

+ 9 2^ 1^ (2fe-l-a)!(a + l)! + 2 ^ ^ (2fe - 1 - a)!(a + 1)! ' 

fc _ m+l Q= V > \ > fc= mil o=m-l V ) \ I 



2fc-l-«/ 



-1) C 



fc=l V 7 fc=l Q=0 

m-1 m-2 n\2k-l-a 



{2k- l-a)!(a + l)! 

1) 



(4.38) 
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Taking into account (4.36), (4.37) and putting (4.32), (4.38) into (4.30) we get 



Co 



E 



\ (2* -1)1 



,h/3 



+ 



-h/3 



2 

m— 1 



m-1 / x \ 

r + E ^^dtr 



fe=i 



+ E ^ r 

fc=l v 



Te h 



m—l 



e h X k X k - e h 

21 m-1 2t-\ 



+ 



2(1 - e h ) 2(1 - e h 



e=i 



m 

TO — 1 



= 1 J 

2£-l 



+ E ^ E <W E E ^+ 

^=1 v ; 1=0 fc=l K i=0 v K y 



m—l 



2^-1 



m-1 x tv 2^-1 



£=1 v ; 7=0 fe=l i=0 K 



A i Q 2^-l-i_ 



7=0 



[m+ll 



2fc-l 



/17 _ 



E E 



(fc/3) 



2fc-l-Q 



(-l) c 



+ 



fc=l a=0 



(2fc - 1 - a)\(a + 1)! 



m—l m— 2 



+ E E 



2fc-l-Q 



i-ir 



(2fc-l-a)!(a + l)! 
P m -2(hP) + 2d- 



+ E E 



(fy3) 2 



-l) c 



i (2fc _ E<W 

1 v ; 7=0 



j' m+l j a=0 v / ; _ |" m+l j a=m 

e -hp e hf3 + e -/i/3 + e l-h/3 + g/i/8-1 _ 4 m -! 



V (^) 2fc 

(2*)! + 



4 E E 



(h/3) 



2fc-l-a 



("1) C 



m—l to — 2 



+ - E E (/ ' ;i)2/ '~ J 



2 ^ ^ (2ife-l-a)!(a + l)! 2 ^ (2fe - 1 - a)l(a + 1)! 

fc=l Q=0 V / V 1 ; /; _ m+1 a= o / v 1 y 



m—l m— 2 



2fc-l-a 



(-1) Q 



+ 2 E E ,2/, - I - )!(n + I j!' 

;,_[m+l a=m.-l V / \ ' / 



(4.39) 



Since the last equation is the identity with respect to (h/3) then from (4.39) equating the coefficients 
of the term (h/3) 2m ~ 2 we obtain that for m > 2 



T = D m (h/3) * f m (h/3) = h. 



(4.40) 



Note that equality (4.40) for some m was proved by directly calculation of the convolution D rn (h(3) * 
f m (h/3) in [23]. 

Taking into account (4.40) for optimal coefficients (4.25) the following formula holds 



m—l 



c f> = h + E { akX i + h ^k~ p ) > = 1,^-1- 



(4.41) 



k=i 
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Keeping in mind (4.40) after some simplifications from (4.39) we get 

m— 1 



oh/3 



C + 



h 



,h 



T + E 



k=l 



a k^h —+ h k- 



'\ k e h -l 



e - 1 
2e 



+ 



fc^ m_1 / \, P h \ N p h \ 1 N 

-<% - + E -A + ^ + 5 E c ^ + M 

fc=l v « « / 7=0 



+ 



/i(l + e 



1 1 c j _ r . v l h>) _ y^ h 2e 
2(1 _ e ft) O0 Z^( 2fc -i)! 2_, 



im-3. 



m— 1 



A fc ^A^ 1 - 



+ E <^ E ^ E ^ E 

^=i v ; j=0 fc=l K i=0 K ■ 



+ E^E^E^E(r^) — 



2^-1 



m— 1 



TV 2t-l 



3=0 
1 2fc-l 



4 E E 



2k-l-ai 



N 



(2fc-l-a)!a! 

/;-[ ">+! ] Q = m— 1 V ' 



7=0 
m-1 2fc-l 



_ e ^(l + e) e-^(l + e) 1 ^ ^ 

4e + 4 + 2 ^ ^ 



2fc-l-Q 



(-1)« 



2 ^ ^ (2fc-l-a)!(a + l)! 

fc_ [ m + 1 ] Q=m— 1 v / \ / 



(4.42) 



As said above (4.42) is the identity with respect to (hf3). From (4.42) equating the corresponding 
coefficients of e hlS , e~ h @ , (h(3) a , a = 0, 2m — 3 one can get the system of linear equations with respect 
to unknowns a,k, (k = l,m — 1), P m _2(/i/3) and d. 



Here we obtain the linear system for unknowns a&, 6^ (k = l,m — 1). 

Equating the coefficients of e ft/3 in the both sides of (4.41) we get the following equation 



m— 1 

E 

fc=i 



(e h -X k )(l- A fc ) K (A fc ^-l)(l-A fc ) 



= 0. 



(4.43) 



Further in (4.42) we consider the terms which consists of (h/3) a , a = m — 1, 2m — 3 and we get the 
equation 



m— 1 

E 



Co 



(21 - 1)! .2. / ^ ( * - ,)! ' ^ (2€ - 1 - ,)!,! £ - ^ (A fe - 1)^ 



+ £ 



E° fc E 



+ 



2£-l 



+ E (w , g :,_ ) - )!]! E'.£ " t ( 1 _tf- 

j=m-l y JJJ k=l i=0 v K; 



= 0. (4.44) 

l - ; _ : . /||/' < ■ ; I -A i. r 1 ' 

J=m-1 

Now from equations (4.3) when a = and (4.4) using identities (4.26), (4.27), (4.28) taking into 
account (4.41) after some simplifications for the coefficients Co and CV we get the following expres- 
sions which are asserted in the theorem: 
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Co 



e h - 1 - h 



m—l 



— i — n \ -> 
h - 1 + 2^ 



fc=l 



' A fc (e fe -e) + Ag(e-l) + Af +1 (l-e fe ) 
fc (e-l)(l-A fc )(e"-A fc ) 



Af +1 (e fe -e) + Af(e-l) + A fc (l-e fe ) 
(e-l)(A fc -l)(A fc e"-l) 



(4.45) 



Cn = 



z h h + l 



m—l 



\N+1 



k=i 



' X k (e - e h+1 ) + Af (e fe+1 - e fe ) + A^+V - e) 
ak (e-l)(l-A fc )(e fc -A fe ) + 



+b k - 



( e - e *»+i) + A£(e ft+1 - e h ) + A fc (e h - e) 
(e-l)(l-A fe )(l-A fc e ft ) 



(4.46) 



From (4.44), using (4.45), grouping the coefficients of same degrees of (h/3) and equating to zero that 
coefficients for unknowns a k and b k we obtain the following m—l linear equations 



m—l 

E 

fc=i 

m—l 

+ E 

fc=i 



22-2 



I 

1=1 v ; ; 

h 2i 



2^ (21-2)1 2-* (A fc -l)»+i 



A fc A'0 2 '- 2 _ \ k (e h -e)+\l(e-l)+\% +1 (l-e h ) 
(e-l)(A fc -l)(A fc -c") 



+ 



3 ,0,-2 2 lz? \ N+i A i 21 - 2 j-x.v.. - 1 >- i 1 - l J - 



E 



(2i— 2)! E (i_A fc )i+i 



E a fe 
fc=i 



A (2Z-D! (x k -iy +1 



A fc (l- e h )+Aff( e -l)+A£ +1 (e*-e) 
(e-l)(A fc -l)(A fc e h -l) 



e fc — 1— ft ft. • 1 [mi 

e^-l 1i J ~ ' I- 2 J' 



2=1 



=0 



m—l 

+ E 

fc=i 



.7 

I 
2=1 



i Af + l A l 2i - 1 

i=0 



^ p^TJT 2- (i-A fc )* +1 



m—l l 
2 J' 



(4.47) 
(4.48) 



Further from (4.3) when a = l,...,m — 2 using identities (4.26), (4.27), (4.28) and the expression 
(4.46) for unknowns a k and we have m — 2 linear equations 



m—l 

E «fc 
fe=i 



'' (l-A fc )'+i A " >- '' S' ^ (1-A fc ) l+1 + " " " 

i=0 2=0 t=0 



m—l 

+ E h 

k=i 



i=0 



(A fc -1) 



- e h l q e 



v 



j-l I 
— 2^,2 V- A fc A 

J ^ (A fc -l) i+1 
2=0 i=0 V ; 



+ 



(e-l)(A fc -l)(A fc -e'') 
\ N k + 1 (e-e"+i )+A fc (e h -e)+A 2 (e^ 1 -e h ) 



+ 



(e-l)(A fc -l)(A fce "-l) 



A 2!(i+l-2)! ft ?TZi— > 3-1, m -2. 



2=1 



(4.49) 

After some simplifications system of equations (4.43), (4.47), (4.48), (4.49) we get the system which 
given in the assertion of the theorem. 

Theorem 4.6 is proved. □ 



From theorem 4.6 when m = 2 we get theorem 4.5. 

For m = 3 and m = 4 from theorem 4.6 we we have the following results. 

Corollary 4.1. The coefficients of optimal quadrature formulas of the form (1.1) with the error 

(3 2) 

functional (1.2) and with equal spaced nodes in the space W 2 ' (0, 1) are expressed by formulas 

l( e h_ e)+A Af (e _l )+Afc(1 _ e fc)\ 




(e-l)(X k -l)(\ k e h -l) 
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k=l v 7 

r _ e h h+l-e» , £ ( n \ k {e-e^)+X^{e'^-e h )+\^+\e h -e) Af +\e~e^)+\j(e^ ~e h )+\ k {e h -e) \ 

^ N ~ eh-1 -l-A,!"* (e_i)(i_A fc )(e*-A fc ) + ° fc (e-l)(l-A fc )(l-A fc e*) I ' 

where a& and 6^ (7c = l,2j are defined by the following system of linear equations 



2 , 2 ,jv+i 

i;=1 fc (A fc -l)(A fc -e*) + ^ ° fe (A fc -l)(A fc e*-l) 2(e"-l) + ( e h_i)2 ' 
2 -jv+i 2 

V nu I- V hu ^ — h ~ 2 I - • 

^ Uk (A fc -l)(A fc -e") + A. °fe(A fc -l)(A fc e"-l) ~ 2(e*-l) + (e*-l) a ' 

(Al^TF + ^E 6fc (A fc -1)* = 12 5 
2 A ^+i 2 X h 

g a fc (A^TF + fc E &fc (Al^TF = T2 ; 
here \k, k = 1,2 are f/ie roois o/ i/ie polynomial 

V 4 (X) = (1 - e 2fc )(l - A) 4 - 2(A(e 2/l + 1) - e^A 2 + l))(/t(l - A) 2 + ^ (1 + 4A + A 2 )) 
which |Afc| < 1. 

Corollary 4.2. TTie coefficients of optimal quadrature formulas of the form (1.1) with the error 
functional (1.2) and with equal spaced nodes in the space W^ 4 ' 3 ^(0, 1) are expressed by formulas 

r _ e h -l-h , V ( n ^(e h ~e)+Xl(e~l)+X^+ 1 (l-e^ £ +\ e h_ e)+X N (e _ 1)+Afc(1 _ efe) \ 

00 ~ e"-l + ^ °* ( e -l)(l-Afc)(e"-A fc ) + (c-l)(A fc -l)(Afce"-l) 



Afc , i A fc _ h—2 j ft 



fc=i v 7 

r _ e h h+l-e h _>4-( n ^k(e-e h+1 )+^ (e h+1 -e h )+X% +1 (e h -e) , a£ + 1 (e-e"+i )+Ag (eW - e *)+Afc (e" -e) 
0JV ~ e' 1 -! fc ^i I * (e-l)(l-A*)(e*-A fc ) + ° k (e-l)(l-A fc )(l-A fc e*) 

where a k and b k (k = 1,3) are defined by the following system of linear equations 

3 3 \JV+i 

, afc (Afc-l)(Afc-e") + °fe(A fc -l)(A fc e"-l) ~ 2(e"-l) (e"-l)* > 
/c=l fc=l 
3 xiv+i 3 

E n , A k l V* h, Afc _ h-2 , h 

; , =1 " fc (Afc-l)(Afc-e") + ^ ° fe (A fc -l)(A fc e"-l) 2(e"-l) + ' 

3 Afe 3 a at+i ^ 

a k (Afc „ fc 1} 2 + g 6 fc (A^Tp = 12 5 

3 a at+i 3 Afc h 

g (Afc fc -l)^ + fc E (Afc-l)^ = 12^ 

3 Afc 3 A^ +2 ft 

E a fc (a„ -i)* + E 6 fc(i-Afc)^ = ~2i; 

fc=l K=l 

k=l k=l 
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here A&, k = 1, 2, 3 are the roots of the polynomial 

p 6 (A) = (1 - e 2h ){\ - A) 6 - 2{\{e 2h + 1) - e h {\ 2 + l))x 



x ( h(l - A) 4 + f (1 - A) 2 (l + 4A + A 2 ) + ^(1 + 26A + 66A 2 + 26A 3 + A 4 ; 



which |Afc| < 1. 



5. Numerical results 

In this section we give the numerical results which confirm the theoretical results obtained in the 
section 4. 

Taking into account (2.9) from (2.21) for the norm of the error functional of optimal quadrature 
formulas in the space we obtain 



(-1)' 



N N 
/3=0 7=0 



sing(/i/3 - h 7 ) f e^-^ - e h ^-W_ _ ^ (/i/3 - fry) 2 *" 1 % 



-2^ 

/3=0 



e h/3 + e ft/3 + e l-hfi + e hP-l _ 4 ^ ( h ffk + (1 _ ^)2fc 



2 

m— 1 

£ 

fc=i 



fc=i 



(2fc- 1)! 



2 • (2k)\ 



+ 



+ 



s 2 - 2e - 1 
2e 



m— 1 



S (2fc + l)! 



(5-1) 

In computations of optimal coefficients we need the roots A&, k = 1, m — 1, |Afc| < 1 of the polynomial 
defined by equality (4.12) 

2m-2 

^ 2m - 2 (A) = P { s m ' 2) ^ S = (! - e 2fe )(l - A) 2m ~ 2 - 2(A(e 2/l + 1) - e h (A 2 + l))x 



s=0 



h(l - A) 



2m-4 



+ 



h 3 (l - A) 2m ' 6 
3! 



E 2 (\) + ... + 



/i 2m - 3 ^ 2m - 4 (A) 



(2m -3)! 



(5.2) 



2k 



where E 2 k(^) = X] e s^ s is the Euler-Frobenius polynomial of degree 2k and for the coefficients of this 



s=0 



polynomial the formula e s = 



j /2k+2\ 



(s + 1 — j) 2k+1 is used which was given by Euler. 



3=0 



For convenience the absolute value of the difference (1.4) of the quadrature formula (1.1) we denote 
by |i?(</?)|. Then by Cauchy-Scwartz inequality we have 



\R( V )\ <Mx)\wt' m - 1] (0,1)|| -\\l(x)\W^ 
We consider the cases m = 1, 2, 3, 4 and N = 10, 50, 100. 

The case m = 1. In the space W^ 1 ' ^, 1) using theorem 4.4 and (5.1), (5.3) for the error of 
optimal quadrature formula (1.1) we have 

N = 10: \R(ip)\ < \\(p(x)\w£'°\o, 1)|| • 0.02886; 
N = 50: \R(<p)\ < \\ip(x)\W^' 0) (0, 1)|| • 0.00577; 
N = 100 : |i%)| < ||^(x)|W 2 (1,0) (0, 1)|| • 0.00289. 



T (rn,m— 1)* 



(0,1) 



(5.3) 
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(2 1) 

The case m = 2. In the space W\ ' (0, 1) using theorem 4.5 and (5.1), (5.2), (5.3) for the error of 
optimal quadrature formula (1.1) we have 

N = 10 : \R(ip)\ < ||¥>(a;)|W' 2 (2 ' 1) (0, 1)|| • 0.000424; 
N = 50: \R(ip)\ < ||p(s)|W 2 (2,1) (0, 1)|| • 0.00001534; 
TV = 100 : \R{ip)\ < ||(^(a;)|W 2 (2 ' 1) (0, 1)|| • 0.37802 x 10~ 5 . 

The case m = 3. In the space W 2 ' (0, 1) using corollary 4.1 and (5.1), (5.2), (5.3) for the error 
of optimal quadrature formula (1.1) we have 

iV = 10: \R(ip)\ < \\cp(x)\W^' 2) (0,1)|| -0.0000108; 
TV = 50 : \R(ip)\ < \\ip{x)\W^ 3 ' 2) (0,1)|| • 0.5643 x 10~ 7 ; 
N = 100 : \R(ip)\ < ||v?(a;)|W 2 (3 ' 2) (0,1)|| -0.6435 x 10~ 8 . 

The case m = 4. In the space W^ 4,3 ^(0, 1) using corollary 4.2 and (5.1), (5.2), (5.3) for the error 
of optimal quadrature formula (1.1) we have 

N = 10 : \R(<p)\ < y{x)\W^' 3) (0, 1)|| • 0.5051 x HT 6 ; 



N = 50 : \R(<p)\ < \\f{x)\W^ 3) (0, 1)|| • 0.3854 x 10" 
7V = 100 : |i?(<^)| < ||vj(x)|W 2 (4 ' 3) (0,1)|| -0.1821 x 10" 



9. 
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